Vortex pinning by a columnar defect in planar superconductors with point disorder 
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We study the effect of a single columnar pin on a (1 + 1) -dimensional array of vortex lines in 
planar type-II superconductors in the presence of point disorder. In large samples, the pinning 
is most effective right at the temperature of the vortex glass transition. In particular, there is a 
pronounced maximum in the number of vortices which are prevented from tilting by the columnar 
defect in a weak transverse magnetic field. Using renormalization group techniques we show that the 
columnar pin is irrelevant at long length scales both above and below the transition, but due to very 
different mechanisms. This behavior differs from the disorder-free case, where the pin is relevant 
in the low-temperature phase. Solutions of the renormalization equations in the different regimes 
allow a discussion of the crossover between the pure and disordered cases. We also compute density 
oscillations around the columnar pin and the response of these oscillations to a weak transverse 
magnetic field. 



I. INTRODUCTION 

A key challenge in the physics of vortex line arrays 
in high-temperature superconductors (HTSC's) is under- 
standing the interplay between vortex interactions and 
various types of pinningi^. The competition between 
thermal fluctuations and pinning can lead to different 
phases such as vortex liquids, Bose and vortex glasses, 
and a more ordered Bragg glas a 3 i 4 i 5 i 6 i 7 . 

Considerable theoretical progress can be made in 
studying two-dimensional superconductors with an in- 
plane magnetic field, where the vortex lines form one 
dimensional arrays. Experimentally this situation can 
be realized using thin platelet superconducting samples 8 . 
The statistical mechanics of such systems is equivalent to 
the physics of interacting bosons in one dimension. The 
low-energy and long wavelength properties can then be 
described within a Luttinger-liquid formalism. For exam- 
ple, one can relate the Bose-glass to vortex-liquid phase 
transition in the presence of disordered columnar defects^ 
to the superfluid-insulator transition in a system of inter- 
acting bosons with quenched disorder 9 . One finds that 
for a given disorder strength the system can be tuned 
across the phase transition by changing the temperature 
(which is proportional to <?, the Luttinger-liquid parame- 
ter). The mapping is different with point disorder, which 
is equivalent to time-dependent point impurities in the 
boson problem. In this case there is a subtle second order 
phase transition between a "supersolid" (with algebraic 
order both in boson and translational order parameter) 
and glassy phase£*i2*ii*i2iiii with decreasing temperature. 

Another important feature of vortex physics in (1+1) 
dimensions is the remarkable response to a single colum- 
nar defect. As argued, originally in the quantum- 
mechanical contexliia and later for vortex arrays, even 
a very weak columnar pinning potential can grow to in- 
finity under renormalization group (RG) transformations 
when g < 1. It was shown that the relevance of a single 
columnar pin at low temperatures leads to a strong sup- 
pression of a vortex tilt induced by a weak transverse 
magnetic field. However, at high temperatures (such 



that g > 1) the pin is less effective, regardless of its 
microscopic strength. Remarkably, the onset of the rele- 
vance (or irrelevance) of a single columnar pin and point 
disorder occurs at the same temperature T* , such that 
9(T*) = 1. 

If both point disorder and a columnar pin are present 
(see Fig.nj then at low temperatures (g < 1) one expects 
a competition between the two: a growing columnar pin 
strength under renormalization leads to stronger corre- 
lated pinning of vortex lines at long length scales. On 
the other hand, the increasing point disorder tends to de- 
stroy the effect of the pin on distant regions. Although 
it was argued that point disorder would always render a 
single columnar defect irrelevant at long wavelengths in 
Ref. ^5|, the precise nature of this competition and the 
different pinning properties above and below the vortex 
glass transition were left unresolved. A detailed study of 
this competition is a primary goal of the present paper. 
We show that in the thermodynamic limit the pinning 
strength is strongest precisely at the transition point T* 
where g = 1. In particular, the number of pinned vortex 
lines is a nonmonotonic function of g and strongly peaked 
at g = 1 . For finite systems the position of the maximum 
is slightly shifted to lower values of g (see Fig.0J). We em- 
phasize that irrelevance of a single columnar pin does not 
imply irrelevanc e of an array of pins. In fact, as it was 
argued in Refs. |3l4j . the vortex glass phase formed by 
point disorder is unstable towards infinitesimal disorder 
in columnar pins, which results in a Bose glass phase. 

The paper is organized as follows. In Sec. [H] we in- 
troduce the model and describe the mapping to a Lut- 
tinger liquid. In Sec. IIIII the renormalization equations 
describing the long wavelength physics are derived. Then 
in Sec. II VI we analyze density oscillations of vortex lines 
near the columnar pin above and below the vortex glass 
transition. At the transition point g = 1, we demon- 
strate that the problem maps onto free Fermions in a 
time-dependent disorder potential. This mapping allows 
computation of density oscillations of vortex lines near 
the columnar pin which we then compare with the renor- 
malization group predictions. In Sec. we discuss the 
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FIG. 1: Schematic view of a planar superconductor with point 
disorder and a columnar pin, represented by a notch cut into 
the sample. The wiggly lines correspond to vortices, which 
alternatively can be thought of as imaginary time world lines 
of bosons in one spatial dimension. 

response of vortex lines to a weak transverse magnetic 
field in the presence of a single columnar pin. Finally 
in Sec. I VII we summarize our results and present conclu- 
sions. Appendix 1X1 estimates the parameter range neces- 
sary to see the effects discussed here, while Appendix IB1 
describes details of our free fcrmion calculation. 



over all possible configurations of vortices weighted by a 
Boltzmann factor proportional to e~~ F °' T . In the limit 
of large sample dimension in a "timelike" direction, Z 
can also be regarded as the zero temperature partition 
function of interacting boson»ii& 

Z = J Du(x,T)e- s , (5) 

with the imaginary time action S = Sq = T$ /T given by 

s ° = ^J dxdr d d - u ) 2 + ^ u f) • ( 6 ) 

For simplicity we here set the Planck's constant Ti = 1. 
Comparison of Eq. (@ with Eq. Q allows us to identify 
the Luttinger-liquid parameter g as 

nT 

The rescalings in Eq. are such that the "sound 

The most relevant contributions to the action from the 
columnar pin S P i n and point disorder S p d reao&iSiiiii 

Spin — V [ drcos[27ru(0,r)], (7) 



II. MODEL 

A one-dimensional array of vortex lines located at po- 
sitions Xj (t) can be described by the density profile (see 

Pig. nj 

n{x,r) =^28[x-x j (T)]. (1) 

3 

Here x and r denote transverse and longitudinal coor- 
dinates with respect to the vortices. It is convenient 
to change variables to the phonon displacement field Uj\ 
Xj(r) = a[j + Uj(r)], where a is the mean distance be- 
tween the vortex lines. In the absence of a columnar pin, 
point disorder, and a transverse magnetic field, the free 
energy of a particular vortex line configuration can then 
be written 

a 2 f 

Fa = — / dxdr [c u (d x u) + c 44 (d T u) 2 ] , (2) 

where en and C44 are the compressional and the tilt mod- 
ulii, respectively. After rescaling x and r 

/ \l/4 /c 4 4\^ 4 
X -> X I — ) , T — > T I — ) , (3) 
\C44/ V C 11/ 

the free energy takes the isotropic form 

^o = ^J dxdr [{d x uf + (d T u) 2 ] , (4) 

with A = a? ! y 'C11C44. The partition function Z describing 
a vortex array at temperature T is a functional integral 



Spd — 2 / dxdTU(x,r) cos[2ttu(x,t) + (3(x,t)] , (8) 

where positive (negative) V and U correspond to repul- 
sive (attractive) potentials. When V < 0, Eq. Q rep- 
resents an attractive columnar pin at the origin. For 
simplicity we take the phase 0(x, t) to be uniformly dis- 
tributed between and 27r and U(x, t) to have a Gaussian 
distribution with the correlator 



U(xi,ti)U(x2 7 t 2 ) = A 5(xi ~ x 2 )5(t 1 - T 2 ), (9) 

where the overbar represents an average over realizations 
of the disorder. The total action S entering Eq. JSJ is 
then the sum of the three contributions 10, J7J, and JSJ: 

S — Sq + Spin + Spd- (10) 

In the following sections we analyze S and various ob- 
servables using the renormalization group. Before pro- 
ceeding with quantitative details we emphasize that al- 
though we focus on the behavior of vortex lines in this 
paper, the action (|10fl can be relevant for many different 
problems such as disordered interfaces^, charge density 
waves which order like smectic liquid crystalsji&iii and 
directed polymer arrays^*. 

III. RENORMALIZATION GROUP FLOW 
EQUATIONS 

Provided one is alert to potential patologies^i, an ef- 
ficient way to analyze disordered problems is to use a 
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replica trick— The noninteracting part of the action © 
then becomes 



So 



-y 



dxdr 



du a dup ^ du 



du 



dr Ot dx dx 



where u a (x,T) is the replicated phonon field and k is 
an off-diagonal coupling which is zero in the bare model 
but generated by the disorder—. It is equivalent to (x- 
and r-dependent) quenched random "chemical potential" 
coupled to the first derivatives of the phonon field 
The replica indices a and /3 run from 1 to n and we 
take the limit n — ► at the end of the calculation^*— 
Equation (|llfl leads to a phonon correlation function in 
momentum space for n — > 0, namely^i 



(u a (k,u})up(k',ci/)) = 



4tt 



k 2 



{g5 a p + «) Sk,k'Saj,u 



(12) 

The other two terms in the action of Eq. i|l(J|) correspond- 
ing to a pinning potential and point disorder, become 



,5' 



(13) 



and 



Spd — 



A J J dxdr cos 2ir[u a (x, r) — up[x, r)] . 

a,/3 



(14) 

To study the statistical physics described by the ac- 
tion H 1 1 p . (|13|l . and (|14|> we employ a momentum shell 
renormalization group scheme 2 ^, where we continuously 
eliminate degrees of freedom depending on frequency and 
momentum within the shell A — <5A < ^uj 2 + k 2 < A. 
Here A ~ l/V a o£o is the ultraviolet cutoff, ao is of the 
order of the lattice spacing, and £o is of the order of super- 
conducting coherence length (single vortex width). The 
resulting renormalization group equations for the run- 
ning coupling constants A (I) and V(l) to leading order 
in A and V are 



dg 

dl ' 
dA 

~dJ 
dn 

~dl 
dV 

~dJ 



o. 



2eA - 2CA 2 , 



= C 2 A 2 , 



V(e- CA) - kV, 



(15) 
(16) 
(17) 
(18) 



where e = 1 — g, I is the flow parameter (A(l) = Ae - '), 
and C is a nonuniversal constant which depends on the 
cutoff A: C a 1/A 2 (in particular, within the shell 
method we find C = 8ng 2 /A 2 ). These equations are sub- 
ject to the initial conditions n(l = 0) = 0, A (I = 0) = A , 
and V(l = 0) = Vb, with A and V being the bare cou- 
plings. Note that the Luttinger-liquid parameter g does 
not change under renormalizatior 



In the absence of point disorder (A(l) = and k(1) = 
0) our results reduce to those obtained by Kane and 
Fisher 1 ^. In this case the columnar pin is relevant for 
g < 1 and irrelevant for g > 1. If V = then our 
equations are equivalent to those first derived by Cardy 
and Ostlund 11 and later extensively explored for different 
problema 5 i 10 i 12 i 13 . Equations (|T3|l - lfT7|) imply that point 
disorder also becomes relevant when g < 1. Contrary to 
bosons interacting with many impurities^ (which is the 
analogue of many columnar defects for a flux problem), 
there is an intermediate fixed point with a finite value 
A* = lim^oo A(Z) = 0(e) which continuously emerges 
from a pure Gaussian fixed point for g < 1. This makes 
the RG approach tractable on both sides of g = 1. Note 
also that when g < 1, n(l) — > oo, suggesting nontrivial 
correlations in this phase. We comment that there are 
some claims questioning the applicability of Eq s. i|15|) - 
(|17|l in the glass phase (g < 1). In Ref. [2{| it was 
argued that a replica symmetric solution becomes un- 
stable for g < 1, resulting in different correlation func- 
tions than predicted by the replica symmetric renormal- 
ization group. However, there is still no evidence show- 
ing that this instability actually occurs. Moreover, nu- 
merical results of Ref. |26j confirm one of the crucial 
predictions of Eqs. (|15fl - l|17|) at g < 1; namely, the un- 
usual behavior of the density-density correlation function 
G{x) oc exp(-ylln 2 |x|). 

In deriving equations (|16|) - (|18fl . we implicitly assumed 
that the cutoff is symmetric in the r and x directions. In 
general this is not true. The anisotropy in the cutoff will 
result in different initial renormalizations of n T and k x . 
However, at large length scales the flows for k t and k x 
look the same and the asymmetry disappears 13 . 

Equations (|15fl - l|18|) contain nonuniversal cutoff de- 
pendent terms. Upon rescaling the disorder potential 
A — > A/C they simplify to 



dg 
dl 
dA 

~dl 

dV 
~dJ 

dn =x2 
dl 



= 0, 

= 2A(e-A), 
= V(e — A) - kV, 
= A 2 . 



(19) 
(20) 
(21) 
(22) 



,ii 



These are the renormalization equations which we will 
exploit throughout the rest of the paper. Note that the 
cutoff A enters Eqs. lfTU|) - (P^ only through the initial 
disorder strength Ao. 

The quantitative predictions of the renormalization 
group equations above are valid if Ao is small compared 
to one. On the other hand if the point disorder is too 
weak, then its effects will be hard to observe in experi- 
ments. In Appendix 1X1 we show that for HTSC supercon- 
ducting films, typical values of point disorder strength 
lie within interval Ao € [0.01,0.1]. This, in turn, im- 
plies experimental relevance of the subsequent analysis 
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of the physics resulting from Eqs. (|19|l - 122|) m the differ- 
ent regimes. 



High temperature phase (g > 1) 

If e is large and negative (|e| 3> Ao) then both point 
disorder and pinning strengths decay exponentially to 
zero as 



A(Z)~ A o e- 2 l e l^0, 
V(l) ~ V e- |e|z -» 0. 



(23) 
(24) 



However, the off-diagonal stiffness k renormalizes to a 
finite nonuniversal value 



.( ^(l-e-^) 



4|e| 



(25) 



As we will see below, the finite value of k(1 = oo) (which 
arises for any e < 0) results in corrections to the power 
law decay of various correlation functions. 



write the asymptotic form of V(l) at large and small I: 



V(l) 



Vbexp(e - A )Z, 
Vexp 



e-A 



I < lo, 
I > / - 

Therefore as I — > oo, V(Z) — > faster than exponentially 
in /. Here, Zo represents a crossover scale 



/o 



2|c| A 



and 



7' = 7 exp 



4 A 4 2 V Ao 



(32) 



(33) 



where Li2(x) is the polylog function 2 ^. In the limit e 3> 
Ao we can use the asymptotic expansion for Li2(x) and 
get 



(34) 



Critical phase (g — 1) 

When e = (i.e. g = 1), point disorder becomes 
marginally irrelevant, and one finds 



A(Z) = 
V(l) = 

«(0 = 



Ao 



1 + 2ZA 
Vb 

(l + 2ZA ) 1 /4 
Ap - A(Z) ; Ap 
2 ^2 



0, (26) 

(An 

e"^ 0, (27) 
(28) 



Note that even though the point disorder A is marginal 
at g — 1, the pinning potential V remains irrelevant (i.e., 
Vil) — > 0) for any nonzero Aq. 



Low-temperature phase (g < 1) 

In the case e > we find the following solutions of the 
flow equations 



A(|) = , 

A + (e-A )e- M 



k(1) = - In 
v ' 2 



l + ^£( e ^_i) 



(29) 

Aq-a(Q 

2 

->e 2 l-*oo. (30) 



Note that k(Z) grows without bound. The explicit math- 
ematical expression for the renormalized columnar pin- 
ning potential is rather complicated. However, one can 



where c is a number of the order of 1. Note that even 
though a term of the order of e 2 appears in Eq. H31JI . 
its presence is justified. Indeed, according to Eq. 122|l . 
at large I we have k oc A* 2 1 ~ e 2 l. It is easy to see 
that higher order corrections in e to the renormalization 
group flow equations (in particular to Eq. l|2U|l ) will result 
in terms of the order of 0(e 3 ) in Eq. (|30|l . 

The parameter Zo defined in Eq. 1)32(1 sets a charac- 
teristic length scale A~ 1 e'°, separating long and short 
length behavior of the pinning potential. As we find be- 
low, it also determines the behavior of various observ- 
ables. Thus, for e > 0, at small I the pinning potential 
first grows under the RG transformations to the value 



V max ~ V(l ) w Vq\J e/ 'A . Then, for larger /, V(l) goes 
to zero faster than exponentially. We comment that for 
g > 1, Zo sets the characteristic scale beyond which A(Z) 
becomes negligibly small and re(Z) stops renormalizing. 

Note that the columnar pin is asymptotically irrele- 
vant in the presence of point disorder for all values of 
g. The mechanisms, which lead to this arc different be- 
low and above the vortex glass transition. Thus in the 
high-temperature phase g > 1 thermal fluctuations are 
responsible for the irrelevance of V at I — > oo. In con- 
trast, in the low-temperature glass phase g < 1 point 
disorder is the cause of the flow of V to zero at large I. 
The distance when the columnar pin starts feeling effects 
of point disorder and becomes irrelevant grows with de- 
creasing Ao. As discussed below, in infinite samples, the 
effect of the columnar pin is strongest (least irrelevant) 
precisely at g = 1. 

In the weak disorder limit one can compute various cor- 
relation functions using the renormalization group analy- 
sis sketched above. In what follows we will discuss several 
quantities of interest. 
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IV. DENSITY OSCILLATIONS AND THE FREE 
FERMION LIMIT 

A. Density oscillations near a columnar pin 

Since the columnar pinning potential is always irrele- 
vant when point disorder is present, it can be treated per- 
turbatively at sufficiently large length scales. The leading 
contribution to the "Friedel oscillations" of the density 
of vortex lines in linear response in Vq is given by 



(Sn(x)) » Vo cos27rnox / f(x, r)dr, (35) 

J — oo 

where the angular brackets represent the thermal aver- 
age and the overbar signifies an average over different 
configurations of point disorder. The quantity 



Sn(x) = n(x) — no 



(36) 



is the deviation of the vortex lines density from the mean 
no = 1/a (see Eq. (QJ). The function f(x,r) is defined 
as 



f(x,r) = (e^K^)-"^. )]). 



(37) 



We note that f(x,r) is proportional to the density 
density correlation function without the columnar pin: 
(Sn(x, T)5n(0 7 0)) oc cos[27mo(x)] f(x, t). Explicitly, one 
finds 



fix T) f» e" 2 /<T dl l9+K(im-Ja(r C -')] 



(38) 



where r = \J x 1 + t 2 A is the distance between the two 
points measured in the units of the original cutoff A. The 
Bessel function Jo (a;) appearing in Eq. I|38(l and in other 
formulas below is nonunivcrsal and depends on the actual 
details of the cutoff procedure. Instead of Jq{x) one can 
use another cutoff function J(x), e.g., a Gaussian, as long 
as it satisfies general requirements J(x — * 0) — ► 1 and 
J(x — > oo) — » 0. If the disorder is absent in Eq. (J3SJ, i.e., 
k = 0, we recover the well-known result for the Luttinger 
liquid f(x,r) oc r~ 2g . If g > 1 then at long distances / 



is given more generally by f(x,T) oc ? 

V = 2[ff + k(oo)}. 



[12|, where 



(39) 



Thus, when point disorder is irrelevant, the exponent of 
the correlation decay becomes non-universal. For g < 1 
the asymptotic expression for / becomesii*i£ f(x, r) oc 

2 i 2 

g— e In x 

Upon using Eqs. (|35|l and (|38|) we find that the behav- 
ior of (Sn(x)) at large distances for g > 1 is 



(Sn(x)} oc Vq 



cos lixn^x 
xT 1 



(40) 



This equation is valid only for x > exp(io), where ^o is 
given by Eq. (|32[1 . For smaller x the exponent r\ changes 
with x. For g < 1 the crossover is now from power law 



decaying correlations of the form (|40() with rj w 2g for 
1^2:^; exp(Z ) to a faster decay 



(Sn{x)) cx V e~ e ln{x) 
in the opposite limit x 3> exp(Zo). 



(41) 



B. Free fermions 

If g = 1 it is well known that using the Jordan- 
Wigner transformation bosons can be exactly mapped to 
spinless free fermions^i. The transformation also holds 
in the presence of a columnar pin and point disorder. 
The columnar pin and the point disorder correspond 
to static and random time-dependent potentials, respec- 
tively. The time-dependent Hamiltonian which describes 
the fermions then reads 



H f {x,T) 



2m dx 2 



+ U(x,t) + V (x), 



(42) 



where U(x, r) is a random potential satisfying 



U(x 1 ,ti)U(x 2 ,t 2 ) = A S( Xl - x 2 )6(t 1 - r 2 ). (43) 

The mass m (corresponding to the tilt modulus in the 
original flux line problem) sets the Fermi velocity Vf = 
kf/m (where kf is the Fermi momentum), which rep- 
resents the sound velocity in the original boson/vortex 
problem. 

If the sample length L T in the timelike direction is 
large, then the partition function of the vortex array is 
proportional to an appropriate matrix element of the cor- 
responding quantum problem^ 



Z = (G\T T e 



J dr J dxli, f (x,t) 



|G>, 



(44) 



where T T is the usual (imaginary) time-ordering symbol. 
This expression is the quantum-mechanical expectation 
value of the evolution operator calculated in its many- 
fcrmion ground state \G), for a given realization of point 
disorder. If the Hamiltonian Tlf is time independent, 
then Eq. (|44|) reduces to the zero-temperature quantum 
partition function. For N noninteracting fermions the 
ground state can be written as a Slater determinant of 
the single particle states. However, because the Hamilto- 
nian Ti.f is time dependent, the states forming the Slater 
determinant will not be the cigenstates of Hf. Instead, 
they will consist of the N largest eigenvalues of the evolu- 
tion operator (see Appendix [5] for further details). Once 
|G) is known one can easily calculate various observables. 
Here we consider the vortex line density 



(n(x)) = -{G\ci(x)c(x)\G), 



(45) 



where c(x) is a fermionic annihilation operator. Suffi- 
ciently far from the boundaries at r = 0, L T , the density 
profile n(x) clearly does not depend on r. 
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Since we are dealing with noninteracting particles, one 
can find the eigenstates of the evolution operator numer- 
ically even in the presence of point disorder. We describe 
details of this calculation in Appendix [BJ Here we just 
mention that we discretize both space and time and write 
the evolution operator as a product of transfer matrices. 
We take a periodic array of M — 201 sites in the space di- 
rection and of size L — 50 in the timelike direction. The 
particle filling factor is taken to be approximately 0.1, 
so that the ground state eigenfunction |G) is the Slater 
determinant of the 21 highest eigenstates. We took the 
odd number of sites to have an exact inversion symmetry 
around the columnar pin in the finite size system and we 
took the odd number of eigenstates to avoid complica- 
tions arising from the double degeneracy of the energy 
spectrum in the absence of point disorder. A colum- 
nar defect of strength Vq = 0.1 is placed in the central 
site, xo = 101. Point disorder is modeled by a uniformly 
distributed uncorrelated random potential on each site 
of the space-time lattice: U(x,t) e [— Uo,Uq], so that 
A = U$/3. For the effective mass in Eq. (I42|l we choose 
m = 5 corresponding to a hopping amplitude J — 0.1 in 
the discretized model (see Appendix 0. For each config- 
uration of disorder we numerically find the ground state 
\G) and the fermion density and finally average over dif- 
ferent realizations of point disorder. In this way we ob- 
tain "Friedel oscillations" of density for different values of 
Aq. In Fig.[21wc plot a calculated density profile of vortex 
lines for Uq — 10~ 3 / 2 corresponding to Ao ~ 3 x 10~ 4 for 
a particular realization of point disorder (top graph) and 
after averaging over about 130 000 disorder realizations 
(bottom graph) . In terms of Ao <|23|) the chosen strength 
of point disorder corresponds to 

3ttA 



Ao 



A 2 



0.3, 



where we used the fact that the cutoff is A w 0.1 for the 
filling factor 0.1. 

Upon fitting the decay of the envelope of oscillations 
(dashed line in Fig. 0) to a power law (see Eq. (|4Tj|) ~) 
for different strengths of point disorder, we extract the 
exponent ij/2. The results are plotted in Fig. [3] Within 
the error bars the dependence of 7/ on Ao is linear as 
predicted by the renoramlization group analysis for g = 1 
(see Eqs. @ and 



C. Boson phase correlations 

Despite the fact that dislocations can not be present 
in vortex arrays (as they are equivalent to magnetic 
monopoles), they can be relevant for other problems, for 
example directed polymer arrays in two dimensions. The 
correlation function, which gives the energy of a dislo- 
cation pair— (or in the quantum-mechanical language 
the boson-boson correlation function^) can be calculated 
similarly to the density-density correlation function 
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FIG. 2: Density profile of vortex line array near a columnar 
pin. The top graph (a) shows the result for a single config- 
uration of point disorder. The bottom graph (b) gives the 
profile after averaging over many disorder realization. The 
dashed line in the bottom graph shows the envelope of the 
decay os oscillations, which is used to extract the exponent r\ 
(see Eq. lUPl). 



where the boson phase angle 4> is conjugate to du/dx. 
Upon integrating out <f> in the standard way, it is easy to 
show that Eq. (|46|l can be rewritten as 



G{x) 



ai[4>(x,r) — (j)(Q,- 



■>i>, 



(46) 



G(x) = ^exp f- J d T u(x',0)dx'jJ. (47) 

The expression above has to be properly regularized to be 
cutoff independent. In the absence of a columnar defect, 
the function G(x) can be straightforwardly calculated as 
in Eq. (J3SJ, yielding 

GWae^r" 111 -^™"'' 1 (48) 

For g > 1 at large x this gives 

G{x) <x x- n/49 \ (49) 

where ij is defined in Eq. I|39|l . This reduces to the 
result in the absence of point disorder when k(1) = 0: 
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Disorder (A Q ) 
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The presence of a transverse magnetic field is man- 
ifested through an additional term to the action in 
Eq. O 10 i^ 



5 n 



nh 
"2ff 



dxdr 8 t u(x,t). 



(51) 



It is easy to see that h is formally a relevant variable 
under renormalization group transformation - in fact its 
recursion relation isifi dh/dl « h. However, when the 
columnar pin is absent the effect of the transverse mag- 
netic held can be eliminated by a shift of the vortex 
displacement field u(x, r) — > u(x,r) — hr. Clearly, the 
contribution to the action (after performing the replica 
trick) coming from point disorder in Eq. I|14(l is invariant 
under this transformation. Thus, the only effect of the 
shift on the action arises from the pinning term 



FIG. 3: Extracted exponent n/2 characterizing decay of den- 
sity oscillations of vortex lines near a columnar pin (see 
Eq. 1401 ') versus the point disorder strength Ao. The dashed 
line represents a linear fit, which agrees with the renormal- 
ization group prediction for g = 1 (see Eqs. 1261 . 1281 . and 



C(ar) oc x x / 2g . For g < 1 as x — > oo we derive 

G(i)«e- t2l "WW. (50) 

Since the columnar pin is always irrelevant in the pres- 
ence of point disorder, it will give only a perturbative 
correction in Vq to Eqs. and (J50). 

Note that the asymptotic behavior of G(x) (Eqs. H49|) 
and JSUJ) is valid only for an isolated pair of bosons (dis- 
locations in the original problem). It can be shownii 
that under the renormalization group bosons are rele- 
vant on both sides of the transition g = 1. Thus, if iso- 
lated bosons (dislocations) are permitted, they destroy 
the supersolid-vortex glass transition. 



V. RESPONSE TO A TRANSVERSE 
MAGNETIC FIELD 

Another way to determine the effect of a single colum- 
nar pin on the flux array is to study the response of 
the vortices to a weak transverse magnetic field h. In 
Refs. and it was argued that in the absence 
of point disorder in the thermodynamic limit an infinite 
number of vortices can be pinned by a single columnar 
defect in the limit h — > 0. While for g > 1 the fraction 
of pinned vortex lines goes to zero with the sample size, 
for g < 1 almost all lines get pinned. In the presence 
of point disorder, because the columnar pin is always ir- 
relevant, the situation is quite different. In fact, as we 
show below, the number of pinned vortices N p (g) has a 
maximum around g = 1. 



Friedel oscillations 



The presence of h modifies the decay of the Friedel 
density oscillations around the pini£. It is easy to see 
that within linear response theory in Vo the expression 
(|35|l transforms to 

/oo 
f(x, t) cos(2nhT)dT. 
-oo 

In particular, for g > 1 the density modulation at large 
distances and weak magnetic fields is readily obtained 
using Eq. (J3SJ: 

(5n(x)) oc y COs27 ^°* (hx)^ K^i(2Trhx), (53) 

where K u (x) is the modified Bessel function of the sec- 
ond kind. One can analyze the asymptotic behavior of 
Eq. at small and large x. Thus, if hx -C 1 (but 
x > 1) 



(8n(x)} cx Vo 



cos 2imox 



0Fr (|) {2nhx) r i- 1 



cos(^)r(^r(^i) 



(54) 

The second term in the expression above gives a negative 
correction to the result obtained in the absence of h. In 
the opposite limit hx 3> 1, density fluctuations decay 
exponentially with x: 



(Sn(x)) oc Vbcos27rnox 



h*- 



-2nha 



(55) 



These findings are consistent with the results of Refs. pjl 
ITt| . where in the pure case when rj — 2g. For g < 1 the 
situation is more involved. Thus, for x <C e l " one can 
still use Eqs. I|54|) and l|55l) with r\ ~ 2g. On the other 
hand if x 3> e'° the asymptotic behavior of (Sn(x)) is 



(5n(x)) cxe- e_ln (x) 



(56) 



for a weak magnetic field, /ii< 1. In the opposite limit 
hx 1 we get Eq. ()55|l with a running exponent rj: 
rj(x) « e 2 In x. 
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B. Number of pinned vortex lines 

A much more striking effect of point disorder on the 
pinning properties of a columnar defect can be observed 
in the number of vortex lines prevented from tilting by 
the defect. This "pinning number" is defined asi^ii 



N p (h) = no 



du(x,r) 

07 



(57) 



where as before no is the mean vortex density. This def- 
inition can be understood as follows: The numerator is 
proportional to the difference between the total imagi- 
nary currents carried by the vortices in the absence and 



in the presence of a columnar pin, respectively. Divid- 
ing this current difference by the average slope h of the 
vortex lines in the absence of the pin gives the effective 
number of vortices which do not participate in the cur- 
rent flow due to the columnar defect. We identify this 
as the pinning number. It is easy to see that the lowest 
order correction to the current due to the pinning poten- 
tial appears in second order perturbation theory. Upon 
using the simple identity 



du(x, t) 



c)t 



= lim • 



3 27ri [u(x,0) — u(x,8t)] 



1 



2ixiii 



(58) 



we can adopt our previous calculational technique to find 



J 



N„ 



n V Q 2 
Anth 



drdr'dx /(0, r — r') [ — j= 



tR(x,t) t'R(x,t') 

j-2 x 2 _j_ ~l 2 



sin [2wh(r - r')] , 



(59) 



where 



R(x,t)= / dl\g + K{l)]e- l J 1 {y/x 2 +T 2 e- 1 ). (60) 



Note that if |e| <C 1 and A <C 1 we can neglect k(1) 
relative to g in the integral to obtain 



R(x,t) 



\J x 1 + r 2 



(61) 



After integrating over x we find 
gVgno 



[ __ h J J drdr' [sgn(r) - sgn(r')] 
x /(0, r-r') sin 27t/i(t-t'). (62) 



The above integral needs to be handled carefully since 
it is sensitive to the order of limits. One way to deal 
with this is to recall that in physical systems the inte- 
gral over the imaginary time is limited by the sample 
size t 6 [—L T /2,L T /2]. Then, in Eq. (1G2H we can make 
substitutions r — > r + £/2, r' — > r — £/2 and use the 
identity 



lim 

L T — >oo 



L T /2 



L T /2 



dr [sgn(r + £/2) - sgn(r - f/2)] = £. 



,( 63 ) 

In this way our final expression for the number of pinned 
vortices becomes 



N„ 



ink 



dUf(0A) sin27r^. 



(64) 



Case I, g > 1. In the high-temperature phase we can use 
the asymptotic behavior /(0,£) oc £,~ r ' to obtain 



V 2 



(65) 



In the absence of disorder, ?/ = 2g and the expression 
above agrees with the one obtained in Ref. |l5| in the 
limit Vb -C h. If the sample size in either the timelike 
direction L T or in the spacelike direction L x is finite, 
then at small magnetic fields l|64|l saturates at 



N p cx V 2 L 3 -\ 



(66) 



where L = minl-L^L,-}. To see this we observe that if 
L T is finite, then the integral over £ in Eq. I|65|) is taken 
within a finite interval |£| < L T . If L x is bounded then 
it is easy to show that the correlation function /(0,£) 
decays exponentially /(0, £) oc e~^/ im for £ > L x and the 
integral in Eq. 165|) is again cutoff at £ s» L x . Note 
that since with point disorder 77 > 2, the fraction of 
pinned vortices always vanishes in the thermodynamic 
limit (L — > 00). 

Case II, g < 1. In the low-temperature vortex glass 
phase, because of the faster decay of correlations with 
distance 



/((U)cxe- 



(67) 



the number of pinned vortices does not diverge as h 
Instead it saturates at 



1/26 



exp 



1 

4? 



(68) 



The first term here comes from relatively short dis- 
tances £ < exp(Zo) and the second one originates from 
£ > exp(Z ), where Z is defined in Eq. J2U). Although 
N p (h) remains finite as h — > for g < 1, note that it be- 
comes very large for Ao <SC e and diverges exponentially 
at e — > at fixed A . 

An important consequence of saturation of N p at zero 
tilt with the system size L for g < 1 is that as L — > 00, 
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FIG. 4: (a) Number of pinned vortices versus the Luttinger- 
liquid parameter g at a fixed point disorder strength (Ao = 
0.2) for different sample sizes L and (b) the position of the 
maximum (<?*) as a function of L. The value of g* approaches 
1 with increasing L. 



the pinning number will have a maximum as a function of 
the Luttinger-liquid parameter at g sa 1. If L is not too 
large (or tilt is finite) then at a given disorder strength 
the maximum of N p (h = 0) will occur at some g* < 1. 
Indeed, if g approaches 1 from below, then N p (h = 0) di- 
verges exponentially (see Eq. JSSJ)). In turn, this implies 
that if L is finite then the system will be insensitive to 
the saturating effect of point disorder if g is very close 
to 1 and the pinning will be quite effective for g < 1 
as in the case without point disorder—. These simple 
considerations agree with numerical evaluation of N p ac- 
cording to Eq. (|64[) . which are plotted in Fig. Q] The 
position of the maximum (<?*) versus the system size is 
also shown in this figure. It clearly approaches g* = 1 as 



L increases. In practice the parameter g can be changed 
varying vortex line density or temperatur e) 15 ! 17 . There- 
fore measuring the dependence of the pinning number N p 
on g should be experimentally feasible. 



VI. CONCLUSIONS 

In this paper we studied the effect of point disorder 
on one dimensional arrays of vortex lines in the presence 
of a single columnar pin. We obtained the renormaliza- 
tion group flow equations for the pinning strength in the 
presence of point disorder and found that the columnar 
pin is irrelevant at the largest length scales (by different 
mechanisms) both for g > 1 and g < 1, which is in con- 
trast to the pure system where it is relevant at g < 1. 
Because the onset of relevance of point disorder also oc- 
curs at g = 1 we were able to quantitatively solve the RG 
equations for g close to 1. In particular we found that if 
disorder strength is weak, then at g < 1 the renormal- 
ized pinning potential V first grows with the increasing 

length scale to V rnax ~ Vo \J (1 — <?)/Ao and then goes to 
zero faster than exponentially. 

Using renormalization group analysis we calculated 
"Friedel oscillations" of the vortex line density near the 
columnar pin. In particular we showed that at g > 1 den- 
sity fluctuations decay as a power law 5n{x) oc \/x n ~ 1 , 
where the exponent r\ is shifted upwards from r\ = 2g due 
to point disorder. If Ao <C <?— 1, then the exponent 77 be- 
comes x independent only at extremely (exponentially) 
large x. At intermediate x one would observe a power law 
decay of 5n(x) with an effective exponent smoothly vary- 
ing with x. In the case g = 1 (which corresponds to dilute 
vortices with a hard core repulsion) we mapped the prob- 
lem to that of free fermions, which we exactly solved by 
direct diagonalization. We found that the decay of 8n(x) 
in this case agrees with the predictions of our renormal- 
ization group analysis. For g < 1 we showed that the 
density oscillations at large distances decay faster than a 
power law Sn(x) oc l/x^ 1 ^ 9 ^ lnx . 

We found that the strongest manifestation of point 
disorder on properties of a single columnar pin is the 
nonmonotonic dependence of the number of vortex lines 
prevented by the pin from tilting (N p ) on g in a weak 
transverse magnetic field. In particular, for g > 1, this 
number diverges at h — > as a power law N p oc h n ~ 3 , 
provided 77 < 3. In finite systems this divergence is cut off 
by the system size and N p saturates N p (h — 0) oc L v ~ 3 , 
where L is the smallest of L x and L T . The only differ- 
ence with the pure case is that rj gets renormalized by 
point disorder ij > 2g. On the other hand if g < 1 then 
arbitrarily weak point disorder leads to saturation of N p 
at h — > 0. As a result in large samples L — > 00, N p (g) 
has a maximum at g* = 1, while if L is finite then the 
maximum occurs at some g* < 1. 
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APPENDIX A: ESTIMATES FOR HIGH 
TEMPERATURE SUPERCONDUCTORS 

In this appendix we estimate effects of point disorder 
discussed above in realistic high temperature supercon- 
ductors. The point disorder correlator A introduced in 
Eq. 10 can be estimated aaLlk 



V4ttA J T 2 w 



3/2 



(Al) 



where $o ~ 2.07 x 10~ 7 Gem 2 is the flux quantum, A 
is the magnetic penetration length, £ is the supercon- 
ducting coherence length, w is the sample thickness, and 
a = 1/no is the one-dimensional density of vortex lines. 
Also j c is the critical current in the presence of point dis- 
order and jo is the pair-breaking current (j c pa jo£, 2 n^ p , 



where 



is the density of point defects). The factor 



of 1/ w in Eq. IjAljl appears after averaging over the sam- 
ple thickness. We are interested at the parameter regime 
where g pa 1. Then using A « A/A 2 (see Sec. lIII|l we get 



\4ttX J T 2 wa \j 



3/2 



(A2) 



Let us take the following estimates valid for typical 
HTSC's: A = 150 nm, T = 50 K, £ = 2 nm, and 
jc/ jo = 0.01. Following Ref. y| we also assume that 
w pa A. Then Eq. I|A2(1 becomes 



A fa 10£/a. 



(A3) 



Quantitative predictions made throughout this paper re- 
quire that point disorder is relatively weak: A <C 1. Al- 
though we expect that arbitrarily large disorder strength 
renormalizes to zero for g > 1 and to small values 
A* pa 1 — g for g < 1 (see Eq. (129(1). it might be nec- 
essary to go to large length scales in order to see this 
effect. From Eq. (|A3|I we see that the condition A -C 1 
is easily satisfied for magnetic fields below H C 2 (we recall 
that in superconducting films with vortices the in-plane 
magnetic field H is related to the interline spacing a and 
the film width w through H pa $o/aw). On the other 
hand, to be observable, A should at least be compara- 
ble to |1 — g\ (see, for example, Eq. (J20J)). If the system 
is tuned within 10% of the vortex glass transition, i.e. 



|1 — g\ = 0.1, then effects of point disorder become rele- 
vant for a < 100 £ pa 200 nm. Such a separation between 
vortex lines is comparable with the magnetic penetration 
length A and thus should be readily accessible at mag- 
netic fields slightly above H c \- 

The other quantity crucial to our analysis is the Lut- 
tinger parameter g. Interesting and nontrivial effects due 
to delicate competition between thermal fluctuations and 
point disorder occur near the vortex glass transition cor- 
responding to g — 1 (see, for example, Fig.^J. In Ref. 
it was argued that g may have a nonmonotonic depen- 
dence on the vortex line density tiq — \/a. In particular, 
at a -C A the Luttinger liquid parameter is very small 
j<C 1 and in the opposite limit a ^> A we have g — > 1. It 
was also argued in Ref. |15| that at least for some param- 
eter values, increasing the vortex line density causes g to 
first increase and then decrease to zero. So one might ex- 
pect that g crosses unity at some intermediate value of a 
of the order of A, i.e., around a few hundred nanometers 
in typical HTSC's. As we showed above, this is precisely 
what we need to achieve point disorder strength of the 
order of a few percent. Any necessary fine-tuning of g 
and A can be achieved slightly varying temperature and 
the in-plane magnetic field, which determines a. 



APPENDIX B: DENSITY OSCILLATIONS OF 
FREE FERMIONS IN THE PRESENCE OF 
POINT DISORDER 

The partition function of free fermions in a poten- 
tial with point disorder (assuming for simplicity periodic 
boundary conditions in the timelike direction) is given by 

Z = Tr(y\T T e~ dT 1 dxHf{x ' T) \y), (Bl) 

where Tr denotes the trace over all possible Slatter deter- 
minants | and Tif is the Hamiltonian (|42() of the non- 
interacting fermions. A slightly different matrix element 
is required for vortices with free boundary conditions and 
finite L T . However, in the limit L T — > oo, this detail is 
irrelevant^. 

If the Hamiltonian were time independent, then in the 
zero temperature limit (where in the quantum language 
temperature is equivalent to 1/L T ) the only contribu- 
tion to the A^-particle partition function comes from the 
Slatter determinant of N eigenstates of Tif with low- 
est eigenenergies (ey). Alternatively these are the states, 
which form the N highest eigenstates of the evolution 
operator 



T = T T c 



- dT f dxH f (x,T) 



(B2) 



Note that, in a time- independent problem, nondegener- 
ate eigenvalues Xj of T become exponentially separated 
from each other with increasing L T : 



X 



3 oc e LT ( £j+1 ~~ Ej \ 



(B3) 
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Therefore, only the largest eigenvalues of the evolution 
operator need to be considered. We are interested here, 
however, in the case where the Hamiltonian explicitly 
depends on time due to point disorder. In this case 
one can not use eigenstates of Hf, because they also de- 
pend on time. However, it is still possible to define the 
spectrum of the evolution operator for a given L T . It 
seems reasonable to assume that even in the presence of 
point disorder the separation between the eigenvalues of 
T grows exponentially with L T m& This statement can be 
checked numerically (see Fig.[5J). Therefore, in the limit 
L T — > oo the partition function i|44|) is dominated the 
"ground state" I*) = \G) of the evolution operator, i.e., 
by the Slatter determinant of the N highest eigenstates 
g l (x) of T: 



\G) = det gi(xj). 



(B4) 



Here i, j = 1, ... ,7V enumerate different single particle 
levels and coordinates. FigureJSJcorresponds to relatively 
strong point disorder and it is necessary to have large 
sample sizes to assure the ground-state dominance. For 
weaker disorder it will suffice to have smaller L T to find 
the exponential separation of \ 3 . 



In order to obtain \G) numerically we discretize the 
system both in space and time directions. In particular, 
instead of the Hamiltonian given by Eq. (|42|) we use 



n f (r): 



M-l 

= £ 

3=0 



[-J( c }+i c j + c ] c j+i) + Uj(r)c]cj +V clc 



(B5) 

Here cj and Cj are fermion particle creation and annihila- 
tion operators, respectively, and M is the number of the 
sites in the spacelike direction. 

We use periodic boundary conditions so that cj + m = 
Cj. The random potential Uj(rk) is taken to be uni- 
formly distributed on each space time point in the in- 
terval [— a, a]: 



Uj(T k )Ui(T p ) = A 5ji8k p , 



(B6) 



where A = er 2 /3. We write the evolution operator as the 
product 




T 



500 



-H f (T k )Sr 



(B7) 



In order to reproduce the continuum limit (|B2|I , we have 
to take the time step 5t small enough to ensure that the 
effects due to [Hf(r), Hf(r + St)] ^ can be neglected. 
However, we do not expect any qualitative changes in the 
results even if this condition is violated. Of course, the 
time step 5t can ultimately be taken to be 1 by an ap- 
propriate scaling of the couplings J, U, and V. In a single 
particle space both the Hamiltonian Hf and the evolu- 
tion operator T are M x M matrices. For a given real- 
ization of disorder, it is straightforward to diagonalize T 
and find its eigenvectors gi(x) using standard numerical 
methods. Then the density profile is given by 



FIG. 5: Seven largest eigenvalues of the evolution operator T 
as a function of the sample size L T . The parameters used in 
this figure are L x = 101, A ~ 0.008, Vo = 0.02, J = 0.1 (see 
Eq. CHOl). For a pure system (Ao = 0), the largest eigenvalue 
corresponding to a bound state is followed by three degenerate 
pairs of smaller ones. Point disorder slightly lifts degeneracy 
between the states, but clearly does not affect exponential 
growth of the separation between the pairs of eigenvalues with 



N 



The index d emphasizes that this is the result for a given 
realization of point disorder. After averaging nd(xj) over 
disorder realizations we obtain (n(x)), which is used in 
Fig. El 
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